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Abstract. Some thermodynamical properties of solids, such as heat capacity and 
magnetic susceptibility, have recently been shown to be linked to the amount of 
ly-j , entanglement in a solid. However this entanglement may appear a mere mathematical 

artifact of the typical symmetrization procedure of many-body wave function in solid 
state physics. Here we show that this entanglement is physical demonstrating the 
Q_i| principles of its extraction from a typical solid state system by scattering two particles 

off the system. Moreover we show how to simulate this process using present-day 
optical lattices technology. This demonstrates not only that entanglement exists in 
' solids but also that it can be used for quantum information processing or for test of 

Bell's inequalities. 



\2 ■ PACS numbers: 03.67.Mn, 03.67.-a, 03.65.Ud 



Entanglement [1] is the property of quantum systems to show correlations that 
cannot be explained by classical mechanics [2]. In the last two decades entanglement 
has been recognized as a resource in quantum information processing, like quantum 
cryptography [3] and quantum computation, [4,5]. Furthermore in recent years the role 
that entanglement plays in the properties of condensed matter systems has received 
an increasing attention. For instance in a typical magnetic solid the electron orbitals 
of different atoms overlap giving rise to an exchange interaction between the electron 
spins which can be described by a Heisenberg Hamiltonian, whose thermal states are 
entangled [6]. 

It is a common belief that entanglement cannot exist on a macroscopic scale 
since decoherence effects would destroy all quantum correlations. However it has been 
predicted that macroscopic entanglement not only is related to critical phenomena [7] 
but also that it can exist in solids in the thermodynamic limit [8] even at high 
temperatures [9]. Indeed there are many works inferring the existence of macroscopic 
entanglement at various temperatures up to room temperature [10-12]. The evidence 
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for this entanglement comes from experiments measuring different thermodynamic 
properties such as magnetic susceptibility and heat capacity. The main aim of this article 
is to show that such entanglement is not an artificial mathematical property but that it 
can be extracted and therefore used for quantum information processing in the same way 
as heat can be extracted and used for work in thermodynamics. In particular we will 
show how such entanglement can be transferred to a pair of particles and subsequently 
used, in principle, for quantum communication or to test the violation Bell's inequalities. 
Furthermore we propose a scheme to demonstrate entanglement extraction with present- 
day technology using optical lattices. In the next two paragraphs we outline the basic 
idea behind our proposal which is then elaborated in the remaining part of the article. 

In the ideal scenario the entanglement between a pair of spins inside a solid can 
clearly be transferred to another pair of probe particles using only local swap operations. 
To this end one sends simultaneously a pair of probe particles toward the entangled 
spin chain in such a way that each probe interacts with a different spin (cf. FigOJ. We 
emphasize that the two probes do not interact with each other nor do they experience 
an interaction mediated by the solid (cf. [13]). The entanglement between the probes 
has been extracted from the spin chain and cannot exist without entanglement in the 
chain. This is a genuine non local process between the two probes like in the case of 
entanglement swapping [14]. We mention that in the continuum limit our procedure can 
be used to extract entanglement from vacuum (some steps in this direction have been 
taken by [15]). In practice, however, the spin chain is in a mixed entangled state and 
the scattering interaction between probes and spins in the chain can only partially swap 
their state. Is entanglement extraction still possible under such realistic conditions? In 
this article we answer affirmatively to this question. We will consider only common 
physical interactions like the Heisenberg or the XY ones. We will show that in general, 
after one scattering event, entanglement swapping is only partial. In this case it is 
interesting to investigate whether by repeating this collision procedure with the same 
probes, one can achieve a better entanglement extraction. 

The most natural way to extract entanglement from entangled electron spins in 
solids would be to scatter pairs of neutrons off the solid. We will present a proposal 
for simulation of this process of entanglement extraction with optical lattices [16] that 
we believe can be implemented at present. Optical lattices are a very useful tool for 
simulating many quantum effects in solids and have been employed in demonstrating 
quantum phase transition [16], and quantum gates [17]. For our purposes we need 
to simulate both the interaction in the solid system as well as the interaction with 
the external probes. Hamiltonians of entangled spin chains or ladders can be realized 
using cold neutral atoms trapped in potential wells generated by counter-propagating 
lasers [18-20]. Bosonic atoms (black in Fig. |2J) are loaded onto the lattice in such a way 
that there is only one atom per lattice site and that the hopping between neighboring 
sites is inhibited. Long-lived atomic states are used to simulate the electronic spin in 
a typical solid. The interaction between atoms can be varied by adjusting the laser 
parameters and by means of electric and magnetic fields. In this way one can produce 
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Figure 1. Schematic drawing of the entanglement extraction process. Two probe 
particles initially prepared in a product state are sent toward a solid whose spins are 
entangled. Each probe interacts with only one spin in the solid. During the collision, 
due to the exchange interaction between probe and spin, their state is swapped. The 
two probes come out from the solid in an entangled state. 



a variety of common spin Hamiltonians. As probe particles one can use different atoms 




Figure 2. Simulation of the entanglement extraction using optical lattices. Cold 
neutral atoms (black) are loaded onto an optical lattice generated by counter- 
propagating lasers in such a way that each potential well is occupied by only one 
atom. These atoms represent the spins degrees of freedom in a typical solid. To the 
first potential a second one is superimposed (wireframe) onto which different atoms 
(white), called markers, are loaded. It is possible to move the markers across the solid 
by moving adiabatically the second potential. This can be efficiently done by varying 
laser parameters of the second lattice. A collision between a marker and a spin can be 
produced by moving the marker into the same potential well of the spin. 



(white in Fig. |2J, called marker qubits [21]. The main idea is to construct first the 
optical lattice, loading the atoms in a register that constitutes the simulated solid and 
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then to superimpose a second optical lattice loading auxiliary atoms onto it as shown 
in Fig. El The auxiliary marker atoms need not to be of the same atomic species as 
the atoms in the solid. The marker atoms can be moved through the solid by varying 
adiabatically laser parameters of the second optical lattice. These operations do not 
affect the register atoms which remain confined in their potential wells. In this way 
collisions between probes and spins in the solid can be simulated by moving a marker 
atom close to a register atom letting their orbitals overlap as shown in Fig. |21 This 
in principle can be performed with high precision and efficiency as in the case of cold 
controlled collisions [17]. In this way one can realize the swap operations needed in our 
protocol. Such an operation can be repeated several, times moving the markers forward 
to collide with other register atoms. 

Let us now introduce the detailed description of the interaction between the probes 
and the spins in the solid. We stress again that all aspects of the following analysis can 
be implemented not only with optical lattices described above but also with any other 
quantum systems capable of coherent manipulation. Two probes are sent to interact 
for a time r each with a single spin of an entangled pair in a solid as shown in Fig. ^ 
We suppose that the interaction Hamiltonian is given by an anisotropic XXZ model 
between probe % = L,R and spin j = 1,2: 

H(X)ij = J (44 + ay y + X44) (1) 

where & x ,y,z are the Pauli matrices, J is the coupling which we assume to be constant 
during the collision. In the case of a time dependent J(t) our analysis still applies 
by making the substitution: Jt — > J* J{t')dt' . The total Hamiltonian is therefore 
Ht = H(\)il + H(X) 2 r. This Hamiltonian comes from the same exchange mechanism 
that gives rise to the interaction between spins in the solid. This Hamiltonian can be 
easily implemented with optical lattices and marker qubits in a probabilistic fashion. 
Here we concentrate on two important limits: A = 1, which is the Heisenberg (exchange) 
interaction; A = 0, which is the XY interaction. We assume that the spin chain is 
initially in the ground state, however our discussion can also be extended to other out- 
of-equilibrium states. We also assume that during the only dynamics is due to the 
interaction with the probes. The reduced density matrix of the entangled pair can be 
written in terms of spin-spin correlation functions by using translational invariance and 
conservation of angular momentum. The full derivation is presented in [22] and the 
result is: 

/ \ + g zz o o o \ 



p 
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2g xx \-g zz 



(2) 



V l + gzz J 

where g zz = (cr*crf)/4 and g xx = (44)/ 4 are spin-spin correlation functions and can 
be numerically estimated for the Heisenberg model. The concurrence of p, assuming 
\9xx I — 9zzi is- 

C = Max[0,-^+4| fe | -2g zz ] (3) 
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Note that the results given below apply to a large class of mixed entangled states 
including, for example, Werner states. Let us first consider what happens after a single 
collision between the probes in a pure product state e.g. 1 00) and the spins when 
A = 1, 0. To illustrate the idea let us suppose that the spins are initially in a maximally 
entangled state \ip~) = 2 1//2 ( 1 01) — 1 10)). After the collision the global state of the four 
particle reads (apart from a global phase factor): 

|¥(f)> = cos 2 Jr \00) lr \r) 12 - i sin 2 Jr \r) LR |00> 12 . (4) 

This state is a superposition of the initial state and of the swapped state. By tracing 
out the two spins one obtains the reduced density matrix p LR of the two probes. We use 
the concurrence [23] to measure the entanglement extracted. For a density matrix of 
two qubits p let us define p = o~ y ® o~ y p*o y ® a y and R = pp. The concurrence is defined 
as C = max{0, Ai — A 2 — A 3 — A 4 } where Aj are the square roots of the eigenvalues of R 
labeled in decreasing order. The concurrence for the two probes is simply C = sin 2 2 Jr 
and is equal to 1 when Jr = 7r/4. Thus no matter how small r is, it is still possible to 
extract some entanglement from the singlet state. 

The extracted entanglement (i.e. the concurrence between the probes), maximized 
over the set of initial product states of the probes, is plotted in Fig. El for A = 1, as 
a function of g xx = g zz , including therefore the case of mixed states of spins in the 
solid. Note that p is entangled for —1/4 < g zz < —1/12 . As it is shown in the 
plot the entanglement between probes decays when g zz decreases. This is confirms the 
intuitive expectation that the more the spins in the solid are entangled, the more such 
entanglement can be transferred to the probes. 

We want to stress that for a translational invariant system the correlation functions 
can never reach the maximum value g xx = g zz = —1/4. Pure maximal entangled 
states can not exist in a homogeneous spin chain. Nevertheless in Fig. El we show the 
concurrence for all possible values of g zz . 

The state of the spins and that of the probes are swapped when Jr = r/4 and 
the entanglement extracted is maximum. In an optical lattice implementation this 
can be achieved in a time comparable to the trapping frequency of the atoms (10/xs 
for 87 Rb, [21]). In the case of solid J is of the order of IK giving approximately 
r ~ 10~ n s [11]. We also note that entanglement oscillates with r, a feature which reflects 
the fidelity of the swapping operation. The probes are always entangled for g zz < —0.16 
no matter how small r is. However, as shown in Fig. |SJ for —0.16 < g zz < —1/12, 
r must be chosen appropriately in order to extract some entanglement. Notice that 
obviously no entanglement can be extracted for g zz > —1/12 since the two spins are not 
entangled. The important point is that as long as some entanglement is shared between 
two spins in the solid (i.e. for g zz < —1/12) our procedure is capable of transferring 
part of it to the probes. Similar results applies also to the XY model A = 0. 

Once r is fixed for a specific experimental setup, the amount of entanglement 
extracted is less than the entanglement in the pair of spins as long as Jr < n/4. Is it 
possible to extract more entanglement by repeating the collision process several times? 
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Figure 3. Extracted entanglement for an Heisenberg chain, measured by the 
concurrence C between the two probes, maximized over the set of initial product states 
of the probes, after one collision as a function of g zz and r. Notice how C increases 
with the correlation g zz and oscillates with the interaction time r. 

This can be done in optical lattices by moving the marker qubits so they interact with 
different identically prepared chains (cf. Fig. EJ). Although this approach has some 
conceptual analogies with the so called homogenization process [24], the two schemes 
are not equivalent since the transformation induced by the interaction Hamiltonian (Q) 
is not a global partial SWAP transformation. Indeed with a local interaction generated 
by the spin Hamiltonian (A = 1) the full transformation of the four spins is: 



U 12L r = PSW lL ®PSW 2R , where (5) 

PSW = e iJt (cos 2 Jt 11 - i sin 2 JtSWAP), and (6) 

SWAP = |00) (00| + |01) (10| + |10) (01| + |11) (11| (7) 

swaps two qubits. 



It is easy to demonstrate that, for A = (XY model), the state LR \^~)i2 is an 
eigenstate of Ht- It is thus a fixed point of the evolution transformation Uulr and the 
state of the probes approaches the state as the number of collisions goes to infinity, 
independently of the input state and of how small r is. This is a counter-intuitive result: 
the state of the probes do not converge to the state of the bath but to an orthogonal 
state. The reason is that XY interaction does not generate a partial SWAP operation 
but something similar to a partial iSWAP operation: 

iSWAP = |00) (00| + % |01) (10| + i |10) (01| + |11) (11| (8) 

that swaps the state of the two qubits but with a i factor. The composition of two of 
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these transformation produces a minus sign that transform in \<P~)- In Fig. 0] we 
show how the concurrence grows increasing the number of collision with the spins. 
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Figure 4. Concurrence between the probes as a function of the number of collisions. 
We supposed the initial state J2J) with 77 = and Jt — 0.2. The data are well described 
by the fit function C(n) = 1 — cxp(— nn) where k = —8.3 10 -3 . 

For A = 1 (Heisenberg model) there is not an eigenstate of the total Hamiltonian 
which is the tensor product of a maximally entangled state of the probes times a 
maximally entangled state of the two spins. 

When the state of the two spins in the solid is mixed, the transformation has still 
a unique fixed point [25,26]. This follows form the Banach-Cacciopoli theorem for 
contractive maps. The transformation we are considering is the tensor product of two 
partial swap transformation which are contractive maps and thus it is itself a contractive 
map. The entanglement of the fixed point depends on g zz ,9xx and r. For example, for 
a Heisenberg interaction, when g zz = —0.2 and Jt = 0.6, the concurrence of the fixed 
point is C ~ 0.28. However, when the two spins in the chain are in a mixed state, 
the entanglement extracted by the probes after repeated collisions is always less than 
that after only one collision for A = 1 (Heisenberg interaction) and "nearly" always, i.e. 
except for mixed states which are very close to a pure state, for A = (XY interaction). 
Therefore repeated collisions do not increase the extracted entanglement. 

Under different hypothesis it is possible by repeated collisions to homogenize the 
state of the probes to the state of the spins in the solid. By modifying the interaction 
between the four spins 1,2, L, R we can build a Hamiltonian that generates a global 
SWAP transformation: 

PSW 1L)2R = e iJt (cos 2 Jt 1 lL ® t 2R - i sin 2 Jt SWAP 1L ® SWAP 2R ) (9) 

This can be found by inspection and reads in the computational basis of Y2LR: 

H 12L r = |0000) (0000| + |0101) (0101| + |1010) (1010| + 1 1111) (1111| 
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+ (|0001) (0100| + |0010) (1000| + |1011) (1110| + |1101) (0111| 
+ |0011) (1100| + |1001) (0110| + /i.e.) (10) 

Notice that the elements of the first two rows can be generated by pairwise interactions 
as in <m>. The elements of the third row instead are four spin processes. With this 
interaction it is proven by the homogeneization process [24] that the probes state will 
approach the state of the spins no matter the initial state. 

Until now we concentrated on extracting bipartite entanglement, however our 
scheme can also be adapted to extract genuine multipartite entanglement. Let us 
consider n spins in a chain, each scattering a spin probe. Let the initial state of the 
chain be a W state: 

\W n ) = (|100 • • • 0> + [010 • • • 0> H h |000 • • • 1» (11) 

and the probes be in state |0) = 1 00 - - -0). W states can be ground states of some 
Hubbard related Hamiltonians [27], critical spin chains [7, 28] and fermionic lattice 
models of high T c superconductivity [9]. The state of the system at time r is the 
analogue for n spins of Eq. (jlj) with \W n ) instead of Thus the reduced density 

matrix for the n probes reads: 

p = cos 2 2Jr |0) (0\ + sin 2 2 Jr \W n ) (W n \ (12) 

Notice that when Jr = 7r/4 the state of the chain is fully swapped onto that of the 
probes as in the case with n — 2. For any value of Jr ^ kir/2, where k is an integer, it 
can be shown that the above state contains genuine multipartite entanglement [28] (i.e. 
it cannot be written as a mixture of biseparable state). 

So far we have assumed that each probe interacts with a single spin in the 
chain. This assumption is, in some cases, unrealistic. For instance it is difficult, with 
present day technology, to address a single spin in a solid. For example, in neutron 
scattering experiments, the spin of the incoming neutron interacts with the total angular 
momentum of a bunch of spins in the solid [29]. It is thus more realistic to analyze a 
model in which two probes interact with several spins at once. Let us consider a chain 
of L spins and let us assume that each neutron, being a wavepacket of a certain width, 
interacts with a different subset of N spins. We will assume that the spin of the probe 
is coupled to the total angular momentum of the N spins, which is equivalent to assume 
that each probe is equally coupled to each of the N spins in the chain. This model is 
equivalent to the ones considered in others contexts [30], and for A = 0, 1 its eigenstates 
are known. Let us consider a Wl state as the initial state of the chain. Though this is not 
the most general entangled state of a spin chain, this model can be solved analytically. 
The concurrence for the two probes for A = is C = sin 2 2 J \^Nt and reaches the 
maximum value ^ for Jr = ty/(4:\/N) (the entanglement for the Heisenberg model 
A = 1 is greater than zero but is always less than that of the XY model). Notice 
that 2N/L is just the fraction of spins interacting with the probes and thus it is the 
maximum bipartite entanglement that can be extracted. We note that these results have 
been obtained for a W state of the chain. This state is symmetric under permutations 
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of any number of spins which means that the entanglement is equally distributed among 

all spins. However for Heisenberg and XY chain models long range correlations are very 

small and the entanglement extracted should decrease with N and with the distance 

between the two beams of neutrons. From the above analysis it emerges that a necessary 

condition to extract entanglement within this realistic model is that segments of the spin 

chains should be entangled. This analysis can be extended the case of non-equilibrium 

states and it would be interesting to see if the entanglement extracted could be larger 

than in the case of the W state. 

We have described a scheme to extract entanglement from spins in a solid by 

scattering probe particles. In this way macroscopic thermal entanglement of solids can 

be converted into a useful resource for quantum information processing. To illustrate 

the scheme we have proposed an optical lattice implementation where all aspects of our 

procedure can be realized with present-day technology. In this context optical lattices 

offer only a simulation of entanglement extraction. This is why we have discussed how 

neutron scattering can be used to achieve entanglement extraction from a real solid as 

well as the physical limitations of this process. We hope that our ideas in the long run 

will lead to an entangling procedure for neutrons in the same way as parametric down 

conversion is for creating entangled photons. 
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